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Abstract 

Several asymptotic expansions of parabolic cylinder functions are dis- 
cussed and error bounds for remainders in the expansions are presented. 
In particular Poincare-type expansions for large values of the argument z 
and uniform expansions for large values of the parameter are considered. 
The expansions are based on those given in |H] and on modifications of 
these expansions given in j^. Computer algebra techniques are used for 
obtaining representations of the bounds and for numerical computations. 
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1 Introduction 

The solutions of the differential equation 

g-(i.^ + .)»^o (1) 

are associated with the parabolic cylinder in harmonic analysis; see [|lO|. The 
solutions are called parabolic cylinder functions and are entire functions of z. 
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Many properties are given in ^ and Q; for applications to physics and many 
more properties see ||] . 

As in 1^ and Chapter 19, we denote two standard solutions of (|l|) by 
U{a, z),V{a, z). Another well-known notation for the parabolic cylinder func- 
tion is Di,{z). The relation between D^{z) and U{a,z) is 

D,{z) = U{-iy-lz). (2) 
Wronskian relations for the solutions U{a, z), U{a, —z), V{a, z) of (|^) are: 

U{a, z)V'{a, z) - U'{a, z)V{a, z) ^ yflf^, (3) 

,dU(a,^z) ,, , , , \/27r , , 

U{a,z) \ > --U'ia,z)U{a,-z)^ / 4 
dz 1 (a + ^) 

which shows that U{a, z) and V{a, z) are independent solutions of (|l]) for all 
values of a. 

Other relations and connection formulae are 



TT 

U{a,z)^ 2 TTi — —r\ [^("' - sin T:aV{a, z)] , 

cos^ nai {a + ^) 



(5) 



T:V{a, z) = T{\ + a) [siuTra U{a, z) + U{a, -z)] , (6) 



27r U{-a, iz) = r(i -t- a) e-'^^^^-i) [/(a, z) -f e*"'*'^-*' [/(a, -z) 



(7) 



C/(a, z) = ie'^^'^C/fa, -z) + , e^"(°'^)t/(-a, zz), (8) 

r(a + i) 

J7(a, z) = -ie-"°f/(a, -z) + p^^i ^ e'^"^°"^^^(-fl, -»^)- (9) 

In III an extensive collection of asymptotic expansions for the parabolic 
cylinder functions as \a\ — > oo has been derived from the differential equation 
(^. The expansions are valid for complex values of the parameters and are 
given in terms of elementary functions and Airy functions. In j^] modified 
expansions are given, which have as extra feature that the expansions are also 
valid when a is fixed and z is large. The coefficients of the modified expansions 
are different from those of Olver's expansions, and they can be generated by 
recursion formulas. 

When Olver published his results, his later work on bounds for remainders 
in asymptotic expansions was not available, and, as he remarked in the 
construction of error bounds for asymptotic expansions of the parabolic cylinder 
functions was an important problem to be considered. In this paper we discus 
error bounds for the remainders of the standard Poincare-type expansions of 
i7(a, z), and of some of the uniform expansions. 
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2 Poincare-Type Expansions 

These expansions are for large z and a fixed. They are given in [Q and derived 
in 0. We have 

[/(a,z)^e-^-%-°-i ^(-ir / |phz|<|7r, (10) 

s=0 '^■^ ^ ' 

^ , 2 , ^ (-a + 7)2^ 



s=0 



By using (||) and (||) the sector of vahdity of (|lj) can be modified, and compound 
expansions are obtained: 



1 ^ (a+^)2s 



s=0 ^ ' 



-e "^-^e^- z"^ ^ , Jtt < phz < 



(12) 



r(a+i) s!(2z2) 



, (a + A)2s 

\/2^ 1,=^ „ 1 (-« + 1)2. 

-^f(^^ ^"^^'^„^!(2;^' -i-<ph-<-4- 

(13) 

With these results we can also obtain compound expansions for y(a, z) for 
other sectors than given in dIT 



2.1 Error bounds of the expansions. 

Bounds for remainders in the Poincare-type expansion follow from where 
results are given for Whittaker functions. The function t/(a, z) is a special case 
of this function. The relation is 

V{a,w) = 2-^^''w-"^Wu,„.{,z), k^-\a, m = i, z^\w^. (14) 
The asymptotic expansion for the Whittaker function reads 

W^fe,,^(z) = zV*^^^+6„(z), «. = (-l)- ^yf , n = 0,l,2,.... 

s=0 

(15) 

We introduce the following quantities. Let 
K 1 
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assuming that a < I. Then remainder e„(z) and its derivative can be bounded 
as follows 



\^n{z)\ < 



2a 



z e ar. 



Vv {t'") exp [2a5Vv {t'^)] , (17) 



V(/) denotes the variational operator (see Q]), which for continuously dif- 
ferentiable functions in a real interval [a, b] is defined by 



VaAf) 



\r{x)\dx. 



(18) 



For a holoniorphic function f{z) in a complex domain the variational operator 
along a smooth arc C parameterized by z{t), a < t < P in which r is the arc 
parameter and z'{t) is continuous and nonvanishing in the closure of {a, P), we 
have ^ 

Vc(/)= / \f[zir)]z'ir)\dr. (19) 

J a 

Along a path V that is a finite chain of smooth arcs (of straight lines, for 
example) V-p can be defined as the sum of the contributions from the arcs. 

In the bounds given in (^) the path V links the point z (with ph z S 
(— Itt, Itt)) to +00, such that on V the condition is fulfilled that di{t + alnt) is 
monotonic. 

We consider the bounds in ([l^ ) for z in the sector [— tt, tt], that is, for w used 
in (|lj) with 5Rw > 0. For other values of w the relations in (||) and ^) can be 
used for computing the function U{a,z). 

In Q simple bounds are given for the variation V- p{t~ ^) appearing in (|l7|), 
for z in certain regions in the complex plane. In Figure 2.1 we show these regions 
in the z— plane and corresponding regions in the w— plane; we only show the 
regions in 5z > 0, > 0, but they should be extended by including the 
conjugated parts. 

z-plane 



V 



Q 



O 



U 



o 



w-plane 


\ ^ 












T 



u 



Figure 2.1. The regions i?i,i?2 and i?4 used for bounding the variation 
Vv (i^") appearing in (p7|). 
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The arc PQ is a circular arc, with radius 2k; Q is the point — -\/3k + «k; S is 
the point in; the arc ST is a circular arc, with radius k. As in ( p^ k = \a\. 

The region Ri is the half-plane Sftz > k; R2 is the region above the curves 
VQ, QS, ST and the conjugated z— values; R4 is the region with \z\ > 2k, < 
K. The corresponding regions in the w— plane follow from z = jw^. 

The following upper bounds for V-p are derived in ||]: 

zeRi, 

X{n)\z\-^, zeR2, (20) 

[x(n) + crw^nl u"|z|"", z e R4, 



XN = 7^ r^, t;= i + Wl-4a2 . (21) 



-1/2 



where 



For Vp {t-^) set n 1 in (|ol). 

For z G i?4, the quantities a, /3 and (5 of (16) should be modified: replace a 
by va and l^l^"'" by u|z|~"'". 

For the parabolic cylinder function U{a, w) these bounds are applicable if 
|w| is large compared with \/[a[. For example, ii w & R4 we need \w\ > 2y^Ja\. 
When using the asymptotic expansion for computations, the restrictions on the 
bounds are not unrealistic because the asymptotic expansion does not make 
sense if w is not large compared with y^jal. 



2.2 Are these domains optimal? 

We have used the same values as in |^ for showing the regions Ri, R2 and 
R4. When verifying Olver's analysis we have found that the boundaries of the 
regions Ri and R2 can be modified. 



2.2.1 Extending region Ri 

For Ri we can use the condition Sftz > max[0, 5R(— a)]. To verify this condition, 
let 

a = u + iv = |a|e*", z = x + iy = re"^. (22) 
In Olver's analysis the path in the 2— plane has to be selected along which 

F{x, y) = ^{z + a\n z) ^ X + \n{x'^ + y'^) - v arctan(?//x) (23) 

is monotonic. Olver chooses an optimal path with the constant argument (j). 
Substituting y — xtaiKp in F{x,y), we find dF/dx ~ {u + x)/x. This is positive 
on the path when SRz > max[0, 3fi(— a)]. One can also use a vertical path x = 
3^0) y 1^ yo upwards, if S (a, a + tt). This follows from 

dF{xo,y) ^ uy - vxq 
dy xl + y^ ' 

Similarly, one can use the vertical path downwards if ^ £ (a — tt, a). In Olver's 
approach vertical paths are not used. 
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Table 1: Ratios p = \€n{z)\/ eViz)- z = re\ r = 10. 



e 








Itt 










§7r 


n ~ 


5 


0.29 


0.30 


0.31 


0.34 


0.13 


0.15 


0.18 


0.25 


0.34 


n = 


10 


0.22 


0.23 


0.24 


0.26 


0.07 


0.09 


0.12 


0.20 


0.37 


n = 


15 


0.18 


0.18 


0.19 


0.21 


0.05 


0.06 


0.08 


0.12 


0.19 



2.2.2 Extending region R-^ 

For i?2 the condition > max[0, 3(— a)] can be used. We verify this by taking 
as Region 2 (as in Olver's approach) the domain with points zq = + iyo from 
which we can draw a half line with the equation xqx + yoy = -\- Vq- This line 
is perpendicular (at z = zq) to the line from the origin to zq. Another equation 
for the line is y = yo — ^^[x ~ xo). On this path we have 

and we see that dF/dx > at z = zq if < ~v < yo or v > 0, yo > 0. With this 
conditions dF/dx > for all x > xq. This explains that we can extend Region 
2 to the domain where Sz > max(0, 5(— a)). 

2.3 Application to the error function 

We applied these bounds for the case a — ^, which corresponds to the error 
function: 

VF_i_i(z) = VttzJ e^erfcV^. (26) 

We computed 

|e«(^)| 

where en{z) is the exact error (see (p^), en^(z) the estimated error (see (|T7|)), 
for several values of n and = phz as given in Table ^ We observe that 
the ratio p is almost always less than i, and that in i?2, where = ^tt with 
J = 4, 5, 6, 7, the estimated error is quite large compared with the real error. 

The upper bound x(»t.) of the variation introduced by Olver is not very sharp 
for certain values of the parameters. In his analysis the estimate of the variation 
along a certain path follows from 

n dr f°° n dr 

|Z -^T^^^p^ " Jo 



where 9 ~ phz and (j) G [—^tt, ^tt] is defined by coscj) — k/\z\, where k — \a\ is 
introduced in (^6|) . The right-hand side in (^8|) is estimated by Olver as follows 

ndT ^ ndT _ X(n) ,^9) 



||z|e»(«-<^') +r|"+i - Jo (|z|2 +T2)i»+i 
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Table 2: Ratios as in Table I, now with variations V-p{t " ) according to (|30|). 
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I- 






8^ 




8^ 






n = 


5 


0.29 


0.30 


0.31 


0.34 


0.38 


0.42 


0.43 


0.41 


0.32 


n = 


10 


0.22 


0.23 


0.24 


0.26 


0.31 


0.35 


0.38 


0.39 


0.33 


n ~ 


15 


0.18 


0.18 


0.19 


0.21 


0.25 


0.28 


0.29 


0.27 


0.17 



The right-hand side of ( p8[ ) can be written as a Gauss hypergeomctric func- 
tion, and we find for the variation (along the same path V) 



Vvit- 



dr 



IzhJo (w2 + 2cos(6i-</.)ii + l)5' 

1 



.^(|:^,;sin^( 



(30) 



This result is also given in j^. The value x(r7,) arises when we replace the argu- 
ment of the F— function by unity. With this new value of V-p we re-computed 
the ratios of Table (|l]), and we give the new ratios in Table (||). We see that 
indeed the ratios become larger in the regions R2 and i?4, except when 9 = tt. 
That is, the estimates of the remainders become more realistic in these regions. 



3 Uniform expansions in terms of elementary 
functions 

We transform the differential equation into a standard form and distinguish 
between the cases that there are no real turning points (as for a > 0), and that 
there are two real turning points. For convenience we consider real parameters. 



3.1 Positive a 



For a > no oscillations occur on the real z— axis. Q gives expansions that 
cover all real z. We consider two different modifications, one for z > and 
another one for z < 0. These modifications are derived in and we take the 
same notation as in this reference. 

The function U{a, z) is a solution of the differential equation (|l|), and w{t) = 
[/(i/x^, /zt-\/2) satisfies 



aw 4,9 , , 

The function W{t) = (t^ + l)iu{yi'^ , i.dV2) is a solution of 

d^W 



de 



w, 



(31) 



(32) 
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where the relation between t and ^ is given by 



^ = Uy/t^ + 1 + i In + 1 



and ip{^) is given by 



(33) 



(34) 



4(i2 + l)3- 

Transformations of this kind are discussed in Chapter 10. The relation in 



3f) follows from the relation 



that Olver used in a Liouville- Green transformation. 
The quantity F in e^^ is a solution of 



(35) 



d^ 



It is convenient to introduce another parameter, t, by writing 

t 



We have 



di 



1 



- 1 



8?2(l + f)2 

and equation (^6|) becomes in terms of r: 



(36) 

(37) 
(38) 



j2 171 jp ^ 

16?2(f+i)2 +r32f(2f2+3f+l)-4^2i +(20f2 + 20f+3)F = 0. (39) 
ciT^ ar 

3.1.1 2; > 

We give an asymptotic expansion of the U— function for a large and positive, 
that holds uniformly for z > 0. We write 



U 



V2/i/i(/i)(i2 + l)4 



where F^{t) satisfies equation (pS) and is expanded in the form 



F,{t)^Y.^-1) 



s=0 



■ Mr) 



where 



h{^l) = 2 



121 12121 



(40) 

(41) 
(42) 



8 



Substituting ( p] ) into ( |39| ) and prescribing 

0o(r) = l, 0.(0) = 0, s>l, (43) 

we find that the coefhcients 0s (t) are polynomials in t of degree 3s, and are 
given by the recursion relation 

0s+i(r) =-4T2(r + l)"— 0s(t)- i / (2Ou2 + 2Ou + 3)0s(w)du. (44) 



For deriving this relation observe that (p3) can be written in the form 



^ dr dr 



1, .idF 



i(20r2 + 20t + 3)F. (45) 



The term h{fi) given in (^2|) follows from (^0|) and from the condition on 
0s (t) given in (||). 

The expansion in ( |4l| ) corresponds with the expansion (11.10) given in In 
our expansion one of the parameters a or t should be large, in Olver's expansion 
a should be large. Both Olver's and our expansions hold uniformly for all real t, 
but we prefer for negative values of i a slightly different expansion that will be 
given in the next subsection. The expansions also hold in unbounded complex 
domains. For details we refer to 

The first few coefficients are 

Mr) -1, 

Mr) =-^(20r2 + 30r + 9), 

Mr) = t^(6160t4 + 18480r3 + 19404t2 + 8028t + 945), (46) 
zoo 

^3 

(f>3 (r) = — (27227200t6 + 122522400r5 + 220540320r4+ 
51o4U 

200166120t3 + 94064328t2 + 20545650r + 1403325), 
For the derivative we have 

f/'(iM^Mt^/2)=-^l±|^l^G.(.), G.(.)^f (-1)^^, (47) 

where the coefficients tps{r) can be obtained by formal differentiating ( ^ ) and 
(HI). It follows that 

Mr) = Mr) + 2T(r + 1)(2t + 1)0s-i(t) + St^t + if ^M^ir) ^ ^^^^ 

dr 
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s = 0,l,2,.... The first few coefficients are 

Mr) =1, 



V'i(t) - -r(28T2+42T+15), 



r2 



V'2(t) = - — (7280T4 + 21840r3 + 23028r2 + 9684r+ 1215), (49) 

V'3 (t) = gY84o (30430400r6 + 136936800r5 + 246708000t4+ 

224494200t3 + 106122312t2 + 23489190r + 1658475). 



3.1.2 z <Q 

For negative z we have 



V / 1(^2 + 2'^'' 



(50) 



where 



S=0 5-0 



Again, these expansions are valid when one or both parameters a and t are 
large. _ _ _ _ 

The functions F^{t),G^{t),P^{t) and Qfi{t) satisfy the following exact rela- 
tion: _ _ _ _ 

F^it)Q^it) + G^it)P^it) = 2. (52) 

The relation in (|^) can be used for obtaining expansions for V{a, z) and its 
derivative. 

3.1.3 Error bounds of the expansions 

We apply Theorem 3.1 of 0, page 366, and write the expansion in (^ij) with a 
remainder. For n = 0, 1, 2, . . . we have 

= E(-l)^^ + RnM- (53) 

5 = ^ 



The remainder RnifJ-, t) can be bounded as follows 

\Rn{fJ;t)\ < exp 



^2 



(54) 
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Tabic 3: Ratios \Rn{n,t)\/ rI^'' ip,t); z = 2ty/E,n^ 3. 



t 


0.0 


1.0 


2.5 


5.0 


10 


25 


50 


a 


= 1 


.21493 


.14455 


.84677 


.94186 


.98360 


.99728 


.99932 


a 


= 5 


.06142 


.43256 


.96494 


.98773 


.99667 


.99945 


.99986 


a 


= 10 


.00343 


.50123 


.98214 


.99382 


.99833 


.99973 


.99993 


a 


= 50 


.04921 


.56597 


.99637 


.99876 


.99967 


.99995 


.99999 


a - 


= 100 


.05601 


.57478 


.99818 


.99938 


.99983 


.99997 


.99999 



where we take into account that we consider positive t and fi (see also Exercise 
3.1 on page 367 of §). 

We need the variation of the coefficients (ps for i > which corresponds with 
^ > and -i < f < 0; cf. @. We have 

rO 



V 







\<l>nir)\dT. 



(55) 



For negative argument the coefficients 4>s oscillate; see Figure |3.1.3 

0.15 




Figure |3.1.3| . Graphs of s = 1,2,3, r e [-1,0]. 

In Table | we show the ratios \Rn{n,t)\/ Rn\fj,,t) (with n = 3), where 
Rn\^i,t) is the right-hand side of (^). We see that the estimates for small 
values of t are much too large. An explanation is that for small t (that is r close 
to —\), the variations in (^ ) are calculated over a larger interval than when t 
is large. 

For z < we consider the remainder in the expansion of (|5l|) and we write 



s=0 



2s 



i?„(/i, t). 



In the present case we have the bound 
\Rn{p,t)\ < exp 



^2" ' 



(56) 



(57) 
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Table 4: Ratios \Rn{fi,t)\/ mf' {^l,t)■, z = -2t^,n = 3. 



t 


0.0 


1.0 


2.5 


5.0 


10 


25 


50 


a 


^ 1 


.29041 


.04469 


.87352 


.76079 


.72513 


.71493 


.71347 


a 


= 5 


.17780 


.02071 


.96996 


.94637 


.93771 


.93509 


.93471 


a 


= 10 


.12433 


.01817 


.98467 


.97279 


.96835 


.96700 


.96680 


a 


= 50 


.07476 


.01644 


.99689 


.99449 


.99359 


.99331 


.99327 


a - 


= 100 


.06829 


.01624 


.99844 


.99724 


.99679 


.99665 


.99663 



where 



^-co#«)= / \<t>ni^)\dr. (58) 



In Table^we again show the ratios \Rn{n,t)\/ Rn\p,t) (with n = 3), where 
in the present case i?!"^' {p, t) is the right-hand side of (p7|). We see that the ratios 
in general are smaller than when z > 0. An explanation is that the variations 
include the contributions from the interval [— 1, — ^J. The parameter r remains 
in the interval [—-1,0], however. 

3.1.4 Upper bounds for the variations of 4>s{t) 

The variations of the coefficients 4>s (?) used in ( p5|) and ( |58|) can be computed by 
numerical quadrature of the integrals, but for real values it is convenient to use 
the zeros of the polynomials 0s(t). For example, 0i(r) has zeros at ti = —0.816 
and t2 = —0.184. Hence, for r G [— i, 0] the variation in ( |55| ) follows from 

The computation of the zeros of (t>'s(T), however, may be not efficient in an 
algorithm. We can avoid this by constructing for the variations upper bounds 
as functions of r. For example, we found for the first few (psir) the following 
simple upper bounds for t S [— ^ , 0] : 



-3t 



V~ {(t)i) < , 

" 4(l-f4.8r^) 

105 ?2 

Vr^(02) < — — , (60) 

32(1 + 18t2) 



-3465? 



3 



128(1 + 52f^) 

The bounds fit at the origin, and are slightly larger at = — i. 



In Figure 3.1.4 we give the graphs of (j)s{T), s = 1 (left), s = 2 (right), and 
s = 3 for r G [— i, 0] . The first graph from the bottom is for (/)s (r) , the second 
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one is for the variation in (^), the third one is for the upper bound given in 
(|60|). The forth one is for the variation in (|5^), and the fifth one is the upper 
bound. Those two are just equal to the second and third plus the variation 
over which equals 0.1692, 0.1602, 0.2415, for s = 1,2,3, respectively. 

Observe that in ( |58| ) we consider re [— ^ , 0] , and that the variation contains 
contributions from [— 1, — i]. 




Figure 3.1.4 . Graphs of (/)s(t), s = 1,2,3, t €E [—^,0], variations and upper 
bounds; see the text for more details. 



3.2 Negative a 

In this case we consider the function [/(— i/i^, /itv^). This fimction satisfies the 
differential equation 



— ^A^^-iK 



(61) 



For expansions in terms of elementary functions three intervals should be dis- 
tinguished: (— oo, —1 — 6], [—1 + 6^1 — 6] and [1 + S, oo) (the turning points 
i = ±1 should be avoided). We only consider the interval [1 + (5, cx)); for the 
other intervals we refer to 19(1 . 



3.2.1 t > 1. 

We write (see 0) 



where F^j^^t) is expanded in the form 



(t2-l)l 



Mr) 



and where h{ij.) is defined in (p2[). 



s=0 



t 



(62) 



(63) 



(64) 
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C = ^Wt^ - 1 - I In , (65) 

and the coefficients (/'s(t) are as in (|l|); see also (||) and (g|). 

The analysis for deriving ( p] ) is similar to that of the case a > 0. The 
function F{t) satisfies the equation (see (39)) 



7dF d 
l-J. — = -4— 



dr dr 
For the function V{a, z) we have 



dr 



- i(20r2 + 20r + 3)F. 



(66) 



V 



i/i2,/i<\/2j 



/iy7r/i(/i)(i2 



— -^P^(i), P,,(t)^5](-1 



where the (ps (r) are the same as in ( ^l| ) . 
For the derivatives we have 

u'[-|^,^^^tV2)^-^hi^,)it^-l)K-^'^G,{t), G,(0^g^ (68) 

and 



(t2 - l)^py« 



Q^(t), Q^(t)^£(-1) 



The coefficients -08 are the same as in (p7|); see also ( ^ ) and 

3.2.2 Error bounds of the expansion 

We write the expansion in (^) in the form 

^MW = E^ + ^:n(M,t)- 



(69) 



(70) 



For the present values of the parameters the remainder i?„(/i,i) can be 
bounded as follows, see (page 366), 



|^«(a*,^)I < exp 



2Voo,c(0i) 



Voo,e(</>n) 



(71) 



In Table || we show the ratios |i?„(/i, t)|/i?,i (/i, i) (with n — 3), where 

R^n\^i, t) is the right-hand side of (|7^). 

From the first few 4)s{t) given in ( |4^ ) and from the recursion relation in 
( PI ) it follows that all coefficients in these polynomials have the sign of (— 1)"- 
Hence, for t > 1, that is, r > 0, the variations in ( [7l| ) can be easily obtained. 
We have for n > 1, using (p^). 



Voo.«(0„) = (-1)^- / ^'n{'y)da = (-1)^</>„(t) = \Mr)\- 



(72) 
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Table 5: Ratios |i?„(/i, t)|/i?l''^ (m, t); z = 2t^,n = 3 



t 


1.5 


2.0 


3.0 


5.0 


10 


20 


50 


a = -l 


.29990 


.57546 


.80676 


.93078 


.98282 


.99572 


.99932 


a = —5 


.69226 


.86898 


.95344 


.98522 


.99652 


.99914 


.99986 


a = -10 


.81624 


.92930 


.97608 


.99256 


.99826 


.99956 


.99994 


a = -50 


.95602 


.98488 


.99510 


.99850 


.99964 


.99992 


.99998 


a = -100 


.97744 


.99236 


.99754 


.99924 


.99982 


.99996 


1.0000 



4 Error bounds by using integrals 

The construction of error bounds of remainders in uniform asymptotic expan- 
sions is available now by Olver's work on differential equations. In this section 
we consider a method for expansion ( pl| ) by using an integral representation of 
the function U{a, z). For this approach it is convenient to concentrate on large 
positive values of z, and to construct an expansion that holds uniformly with 
respect to a € [0, oo). This expansion reduces to the Poincare-type expansion in 
10|) when a is fixed after expanding the quantities in the expansion that depend 



on A = a/z for small values of this parameter. In fact, by writing /i 



2A 



and t = 1/(2VA) the same can be done in the expansion given by (|4^) and (^. 

4.1 An integration by parts procedure 

We summarize from B and start with the integral representation 



1 ^2 
"1^ 



U{a, z) 



r(a + i) 7o 
which we write in the form 



— ■t;W—ZW dW 



a > — T 



U{a,z) 



nil 1^2 

r(a + i) 



(73) 



(74) 



where 



[w) = w + ^w^ — Alnw, A = — = — r, 



(75) 



4t2 

where t is used earlier in the notation [/(^/i^, /ii\/2). The positive saddle point 
Wo of (f){w) is 

' (76) 



Wo = 5 



V1 + 4A- 1 . 

A standard form of ( [74| ) is obtained by using the transformation 

(j){w)=s-Xliis + A, (77) 
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where A does not depend on s or w; we prescribe that w = should correspond 
with s = 0, and w = wq with s = A, the saddle point in the s— plane. This gives 



A ^ ^Wn + wq — WnwQ — X + \ ln X, 



U{a, z) 



where 



/(s) = (l+4A) 



2^0 

z e 
(l+4A)3r(a + i) 7o 

dw 



\ I s 



ds 



= (1 + 4A)' 



s-X 



s w'^ + w — X' 



(78) 
(79) 

(80) 



By normalizing with the quantity (1 + 4A) we obtain /(A) = 1, as can be 
verified from ( |80| ) and a limiting process (using I'Hopital's rule). 

For A ^ the saddle point wq tends to zero, and the mapping becomes 



It is not difficult to verify that for A = we have 



'1 + 



V 2(1 + 2s) ■ 

If A 7^ the transformation (|77| ) can also be written 

Au) = u;osei(^"''+"'-^-5"'o-'"''+^\ 

in which form no logarithms occur. 

We introduce a sequence of functions {fk} with /o(s) = /(s) and 



/fc+i(s) ^ Vs-r 
ds 



- fk{s)^fk{X) 

* 7~x 



k = 0,1,2, 



(81) 



(82) 



(83) 



(84) 



The expansion in ( ^ ) can be obtained by using an integration by parts 
procedure. Consider the integral 



We have (with A = a/z^) 



r(a + i)7o 



Fa{z) 



-2a-l 



/(A) 



V(A) 



7(A) 



r(a + i) Jo 



z^r{a+^) Jo 



z^r{a+^) Jo 



2 fJ H 

s-e— /(s)^, (85) 



2 S 

s-e- •[/(,) - /(A)] ^ 



s~X 

\ s 
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where /i is given in (p3) with /q = /. Repeating this procedure we obtain 



U{a,z) 



■E 



fkiX) 



z (1 + 4A)3 fc=o 



2k 



(86) 



The factors in front of the series in (Bq) and (^l|) are the same. This can be 
verified by using a — and z = fj,v2t. Also, the series correspond termwise 
with each other, the relation between the coefficients being 



</.fe(f) = (-l)'=(2A)^-A.(A), T 
For example, we have 



1 



V4A+ 1 



- 1 



/o(A) = l, /i(A)--^pl^(20f2 + 30f+9). 



24(f + 1) 



(87) 



(88) 



We write (pq) with a remainder: 



U (a, z) = 



where 



z ^(1 + 4A)3 
i?„(a, z) : 



E 

.A;=0 



r(a + Jo 



^Rn{a,z) 



(89) 



(90) 



4.2 Bounding the remainder 



From ^ and (fg) we infer that /(s) = ©(s^i/^) as s oo (see also (||) for a 
simple verification when A = 0). So, /(s) is bounded on [0, oo). Further, we can 
prove (which will not be done here) from (84) that /„(s) = 0(s^^/^), n > 1. 
Hence, the functions /,i(s) can be bounded |/,i(s)| < M„(A), for s > 0. Using 
such a bound in ( pO| ) will indeed give an upper bound for \Rn{a,z)\, but this 
bound may be not realistic. 

A much better bound will be obtained by estimating |/n(s)| accurately in a 
small interval around s = A, and accepting a rough estimate for other s— values. 
This can be achieved by using a "weight function" w„(s. A), and by writing 



\fn{s)\ < [|/„(A)|+M„(A)] w^{s,X), 
where, for example, we take 

u;„(s,A)= [(s/A)-^e^~^]"". 



(91) 
(92) 



We have Wn(A, A) = 1, that is, for the s— value where the dominant part s°'e ^ * 
of the integrand in (pW assumes its maximal value when a and z are large. We 
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try to find M„(A) > and cr„ > such that (|9 
obtain the bound 



holds for all s > 0. Then we 



\Rnia,z)\ < [|/„(A)| +M„(A)] Snia,z), 



where 



(7n\^'^"^''^^ r(a + 5 - Act,, 



r(a+i) 



For Sn{a, z) we need the conditions 



z'^ > a„ 



a + I > Act„ . 



(93) 

(94) 

(95) 



The quantity S'„(a, z) is close to unity when a + z is large; because /„ is bounded 
the value of ct„ will be small. Numerical calculations show that for ct„ = 1 and 
z > 3, a > 1 the maximal value of S'„(a, z) is smaller than 1.062 . 

In Figure we show the graph of (1 + 5A) [|/„(A)| + A/„(A)] (the factor 
(1 + 5A) is chosen because of scaling) when we take cti = 1 in (|9l|). We have 
Mi(0) = 0, /i(0) = -3/8 and /i(A) - 1/(48A) for large A. We also draw the 
graph of 

\Ri{a, z)\ 



PI (A) 



(96) 



[|/i(A)|+Mi(A)] 5i(a,z)' 

the ratio of the exact error and the estimated error. We see a sharp dip at 
A = 8.3176..., which is a zero of /i(A); for this value of A the asymptotic 
approximation improves, as expected. For large A the quantity pi tends to 1. 

We computed Mi in [|/„(A)| + M„(A)] as the inf |/i(s)/w;i(s, A)|, s > 
and this gives a continuous function Mi (A), which need not to be smooth. 
For example there is a noticeable non-smooth behavior near A = 11, because 
inf |/i(s)/?i;i(s, A)| occurs for A < 11 in a different s— domain than for A > 11. 



1.0- 



0.75 



0.50 



0.25- 




X 40 



Figure ^ The graph of (1 + 5A)[|/i(A)| + A/i(A)] with Mi(A) introduced in 
(|l]) (with CTi ^ 1), and pi(A) from (||). 
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4.3 Bounding the remainder by using Cauchy-type inte- 
grals 

Computing the functions /n(s) by using formula is quite difficult, also 
when we use computer algebra. The representations contain derivatives and 
removable singularities at s = A. In particular the poles are very inconvenient 
when computing the functions /n(s) near s = A. 

It is possible, however, to represent fn{s) as a Cauchy-type integral. The 
mapping in (^) is singular at w = w- ~ — I — wq, the negative saddle point of 
4>{w) defined in (|75|). The corresponding s— value s_ is negative. If A = then 
W- = —1 and the corresponding s— value is s= — i. See (|82|), where indeed /(s) 
shows a singularity at this point. For large values of A we have the estimate 
(see @, formula (4.45)) 



s(u'_) 



-A 



0.2785 



0.4356 



VA 



(97) 



For constructing the Cauchy-type integrals we use the property that the 
functions fn{s) are analytic functions in a domain V in the half plane 3?s > s_; 
in particular, V contains a neighborhood of the positive real axis. 

As in [Hi, we start from 



If 1 

fnis) = — / (3o(ct, a, s)/„(cr) da, Qa{(T, A, s) = , 



(98) 



where s € V and C is a contour in V around the point cr = s, we obtain by 
using the recursion (B4h 



(99) 



where C is a contour in V around the points a = X and ct = s, and 



1 



d 

Qo + 2cr—Qo 
oa 



2{a - A) 

Continuing this we obtain for n 0, 1, 2, 



2((7- A)(cr-s)2 



fn{s) 



1 

2Td 



(3„(ct, A, s)/(cr) da, 



(100) 



(101) 



where C is a contour in V around the points a — X and a — s. The rational 
functions Qn follow from the recursion relation 



Qn — 

For example, we have 

Q2ia,X,s) = 



1 



2{a - A) 



d 

)„-l -I- 2a TT-Qn-l 

oa 



3a^ - Ao-2 - as^ - Xs^ + 6a^s - 6Xas 
4(cr- A)3(a-s)3 ' 



(102) 



(103) 



19 



The coefficients /fc(A) in ( p6| ) follow from (101) by substituting s = X. 

In order to obtain a bound for /^(s) with representation in (101), we select 
a special contour. We take for C the vertical hue — — ctq, where (Tq > 0. 
First we write the quantities Qn in a special form: 

Qi(.,A,.) = ^^-ii±^^^--^ + l (104) 
2p^q q 2pq 

where 

p = a — s, q = (J — X. (105) 

Similarly, 



8A^g + 16Xq^ + 8g3 + 4pA^ - 4pq^ - 2Xp^ -~ p^q 
Q2{(J, X, s) = — ^-^ , (106) 



which also can be written as a sum of partial fractions in which a occurs only 
in the denominator. 

For (T G C, we write a = — (Tq + «t, t G M. For s > and A > 0, we have 



IpI > kl - y^g+T^, \q\ > \X + ao + iT\ > ^al + T^. (107) 
This gives for Qi the bound 

((To+t)' 2(o-o+T^) 

We also need a bound of /(s) for s G C (we know that this function is bounded 
on C, see beginning of Section Let 



M(ao,A) = max|/(s)|. (109) 



This gives for the remainder defined in (|90|) the upper bound 

4A + Sttcto 

Similarly, 



Ri{a,z)\<M{ao,X) ^ , (110) 



,^ , ^ ,, SA^ llA 13 

A, .)| < 2(^2+^2)5/2 + 4(ag+r2)2 + sRT^' ^^^^^ 



9A^ + UttActo + 26cr^ 



|i?2(a,^)| <M(ao,A) ^ ^, (112) 



15A3 eiA^ 43A 81 

iy3(<^>A,s)| < (^2 + ^2)7/2 + 2{al + r2)3 + 2(^x2 + r^f/^ + 8(^2 + r^f ' ^^^^^ 
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|i?3(a,z)|<Af(ao,A) ^ ^ ^. (114) 



reSA-^ + 5497rAVo + 1376Ac7g + 2437rcrg 

The singularity s_ estimated in ( |97| ) is of order ©(A), and it follows that ctq 
can also be taken of order 0{X). When we choose cto in this way, the quantities 
Qm for this choice of C, are of order ©(l/A""*"^) as A ^ oo, and the remainders 
are of order 0(1/A"). The coefficients /^(A) also are of order ©(l/A"). This 
follows from ( |87| ) and properties of 4>k{^)- 

In the above method the estimates of higher Q„ are quite easy to obtain; 
only the estimate of |/(s)| is needed for obtaining bounds of Rn{z,a). From 
numerical verifications we conclude that |/(s)| is maximal on C at s = —(To, 
that is, we can take M{ao, A) = |/(— cto)|- Also we conclude that the bounds for 
the remainders obtained with this method are less realistic than those obtained 
by the other methods, unless A is becomes large. 



5 Numerical aspects 

For computing the numerical upper bounds of the remainders in the expansion 
we used computer algebra for manipulating the formulas. This became already 
quite complicated, although we didn't consider complex parameters so far. It 
was necessary to develop new algorithms for the parabolic cylinder functions, 
because for certain cases we needed accurate values of the asymptotic expan- 
sions and the function U{a,x) for computing values of the remainders in the 
expansion. 

For example, for computing Rn{fJ,,t) of ( p3| ) for n = 3, a = 100 and t — 50 
(which corresponds with x = 1000), the values of F^(i) and the asymptotic 
series in ( p3| ) (with n = 3) arc 

0.99999962523819834461, 0.99999962523819834799, 

respectively, with 17 corresponding digits. We have computed these values with 
Digits = 30 (in Maple), and developed algorithms with adjustable precision, 
based on quadrature methods. We will present these Maple codes for com- 
puting the function U{a,x) for real parameters on the web si te of our project 
( ^ittp : / / tur ing . wins . uva . nl^thk/ specf un/ compalg . html ) . Fortran versions 



of these codes will be developed in a different project. 
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